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Abs tract 

The effect of the Coulomb interaction in superconductivity 
is investigated by means of the self-consistent method employed 
in the boson formulation of superconductivity. The charge and 
current density of both quasi fermion and boson parts are 
calculated. It is found that the charge and current density in 
the ground state are dominated by the boson part. Also, our 
study seems to indicate that the Goldstone theorem is still 
valid, even in the presence of long range forces in the sense 
that a boson condensation is permitted. This massless particle 
recovers the gauge invariance and is responsible for the 
dynamical rearrangement of symmetry. 
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1) Introduction 

It is well known in the theory of superconductivity that 
the presence of collective modes is essential in order to have 
a theory which is implicitly gauge invariant. In the original 
BCS formulation the derivation of the Meissner effect was not 
gauge invariant; this feature was soon understood as due to an 
unsatisfactory level of the perturbation methods used. Several 
calculations^ performed in approximations of higher order than 
a simple Hartree-Fock approximation showed how at this level 
the gauge invariance of the theory is restored. At the same 
time, in higher order approximations, collective states which 
have been interpreted as bound states of quasi-electrons with 
opposite momentum and spin appear. These collective modes, 
which behave as boson fields, restore the gauge invariance and 
the conservation of the current. In this sense the presence 
of these modes is essential in the construction of a theory of 
superconductivity which is naturally gauge invariant at every 
stage of the calculation. 

This situation has motivated a different approach, the 
boson method in superconductivity, which has been developed in 
a series of articles * 2 ' 3 ' 4 * . The attitude of this approach 
is based on the fact that the role played by the boson field 
should be investigated at every level of the theory. In the 
course of our study we have found that the inclusion of the 
boson field besides recovering the gauge invariance/ enables us 
to formulate the theory in a systematic and unified form which 
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offers several practical advantages when applied to particular 

(5) 

problems, such as the appearance of flux lines in the mxxed 

( 6 ) 

state of type-II superconductors and Josephson effects 
Indeed the boson method turns out to be of practical applicability 
in those regimes where the Gor'kov equations cannot be reduced 
in the form of differential equations (Ginzburg-Landau equations 
and their generalizations)^. In preceding articles ^ 3 ' ^ the 
charge and current density has been assumed to be principally 
described by the boson field; the quasi-fermion contribution 
has been considered negligible. One of the results of the 
present article shows that this assumption is correct. 

A study of these collective modes is also important from 
a more fundamental point of view in the phenomenon known as 
spontaneous breakdown of symmetries. There exists in nature a 
variety of many-particle systems which share in common the 
feature that the observed behavior does not present all the 
symmetry properties possessed by their underlying Hamiltonians . 
In general these systems have a ground state which presents a 
systematic structure and is asymmetric with respect to the 
symmetry transformation. In the case of superfluid systems 
(superfluid Helium and superconductors) the Hamiltonian is 
invariant under gauge transformations, but this symmetry is 
phenomenologically violated. It is a fact that a collective 
oscillation whose energy spectrum vanishes in the limit of 
long wave lengths is always associated with the presence of an 
asymmetry in the ground state. This situation has been proved 
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in complete generality in the well known Goldstone theorem 
It is also known that this theorem does not provide a complete 
description of the situation since the fundamental question of 
what happened to the original symmetry remains unanswered. 
Since the invariance of the theory must be kept at every stage, 
it may happen that owing to the highly non linear character of 
the interactions the symmetry takes a different form when the 
theory is expressed in terms of quasi-particles (incoming fields 
in the terminology of quantum field theory) . The conservation 
laws will still hold at the level of the physical fields (observed 
particles) , only they will show different forms according to the 
new shape taken by the symmetry. This phenomenon has been 

(9) 

denoted by the name of dynamical rearrangement of symmetries 
In the process of rearrangement of symmetry the Goldstone 
particle plays a principal role; in the course of numerous 
calculations performed on different models it has always been 
found that the Goldstone particles rearrange the symmetry * 10 * . 

It is now a fact that the Goldstone theorem does not appear 
consistent with the experimental situation in superconductivity; 
this circumstance has led one to expect that the Goldstone theorem 
may break down when long range interaction of the Coulomb 
type are present. In particular let us consider the case of 
superfluid liquids and superconductors; below a certain 
critical temperature both systems exhibit the phenomenon of 
condensation, which implies an asymmetric ground state and 
spontaneous breakdown of symmetry. It is phenomenologically 
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observed that in superfluid liquids there are collective modes 
whose energy goes to zero in the long wave length limit, while 
in superconducting systems such modes possess a finite energy gap. 
The only sensible difference which exists between these two 
systems is the presence of a Coulomb interaction among the 
particles. The disappearance of the Goldstone particle and the 
presence of long range interactions were theoretically related 
by Anderson . He found that a neutral system of fermions with 
an attractive short range interaction exhibits a collective mode 
whose energy vanishes in the limit of long wave lengths; when 
a Coulomb potential is introduced, these collective excitations 
acquire a finite mass which coincides with the value of plasmons 
which are usually observed in many-electron systems. 

Other examples of a failure of the Goldstone theorem have 
been studied in relativistic models containing gauge f ields ' 

However all these theoretical investigations do not provide 
an answer to the fundamental question: how is the symmetry 
rearranged without the Goldstone particle ? 

In the present article we consider a superconducting system 
taking into account the effects of the Coulomb potential. We 
are particularly interested in studying how the presence of 
long range correlations influences the collective modes. This 
problem is of interest at two levels. On one hand it is 
important in the framework of the boson method since the central 
point (4) of this method is the introduction of some invariant 
transformations expressed solely in terms of the boson field. 
I 
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By applying these transformations we can describe space and 
time dependent situations^ 5 ' 6 *. On the other hand this problem 
is of general interest in the context of the spontaneous break- 
down of symmetries. In particular we wish to understand the 
effect of long range forces on the Goldstone theorem. Does 
the Goldstone particle completely disappear from the theory? 
How is the symmetry rearranged at the level of the physical fields? 
These questions have been left unanswered in previous treatments, 
but in our opinion they play a basic role. In this article 
we examine these problems and our study gives a very clear answer 
to these questions. 

(13) 

By applying a self-consistent calculation in the pair 

approximation to the well-known Hamiltonian for electrons in 
superconductors we express the operator form of most of the 
observable quantities in terms of quasi -electron and boson 
fields. We prefer to express the theory in terms of these 
operator forms instead of the usual approach based on the Green 
functions because these operator expressions permit us to 
perform those invariant operator transformations (boson trans- 
formations) which represent the cornerstone of the boson method. 
Furthermore these operator forms permit us to investigate the 
problem of the rearrangement of the symmetry. 

In section 2 we derive the equations which fix the 
coefficients of the dynamical map. These equations are solved 
in sections 3 and 4 where explicit expressions of the current 
and charge density in terms of the physical fields are derived. 
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Finally in section 5 we discuss the results obtained in 
connection with the Goldstone theorem and the dynamical re- 
arrangement of the symmetry. We obtain an equation for the 
boson field which differs in a sensible way from other equations 
obtained in different approaches (1) in the sense that it permits 
not only excitations at plasma level but also condensation in 
the ground state. This condensation carries the gauge trans- 
formation and is responsible for the invariance of the theory 

at the level of the physical fields. The presence of this 

(5 6 ) 

condensation is also important in practical applications 
since it permits us to treat time- and space-dependent phenomena 
in superconductivity avoiding the different problem of solving 
the Gor'kov equations. 

In closing this section we would like to mention that 
the problem of Coulomb effects has been already investigated 
in Ref. 3; however, the treatment is not completely satisfactory 
since the effect of the quasi-f ermion contribution to the charge 
density was completely neglected. In this article we apply 
the self-consistent method to the BCS Hamiltonian including 
the Coulomb potential; we find that the expressions of the 
quasi-fermion contributions to the charge and current density 
are strongly suppressed by the Coulomb potential with respect 
to the neutral case. As a consequence of this result the boson 
transformation, introduced in Ref. 4, turns out to coincide with 
the original gauge transformation for all the frequencies of 
physical interest. 
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2) Dynamical Map 

A superconducting system in the absence of magnetic fields 
is described by the well-known Hamiltonian: 

H = H + H, + H (2.1) 
o 1 c 



H = Z ( d 3 xi|. + (x)e(V 2 )iJ> (x) (2.2) 
o s J s r s 

H 1 = A | d 3 x^| (x) f[ (x) ^ (x) (x) (2.3) 

H f d 3 x f d 3 y P J^M (2.4) 



where 



e(V 2 ) = --^(V 2 +k 2 ) (2.5) 

p(x) = ^(xJ^U) + i^|(x)^(x) (2.6) 

and s = + , + . Following the methods of Refs. 13 and 2 we shall 

assume that the equilibrium state of a superconductor can be 

described in terms of a set of free fields {* ,B} ; corresponding 

to the quasi-electron fields; B corresponds to a collective 

field. It should be noted that the existence and the properties 

of these fields are not assumed; they are the results of the 

dynamics of the system described by the Hamiltonian (2.1) and 

will be determined in the course of a self-consistent calculation. 

Introducing the corresponding set of creation and annihilation 



operators {ol^ + (t) ,a k+ ^ (t) } and {bJ (t) ,B k ( t) } we shall con- 
struct a Fock space as a realization of the canonical commuta- 
tion relations: 

I a ks'V J + = 5 rs 5( ^ (2 * 7) 
lB k/ B^]_ = S(k-p) (2.8) 

We shall not introduce any other Fock space; every other operator 
present in the theory should be expressed in terms of the free 
operators. Thus, introducing the annihilation and creation operators 

a kt,*' a k + ,i ° f the fieids 

we shall express them as a linear combination of normal pro- 
ducts of the free operators: 

a ki (t) = cos8 k a k+ (t) - sin9 k a! k+ (t) (2.10) 

here a and ot + are intermediate operators expressed in terms of the 
physical fields. Since in the following we are going to use only 
bylinear products of such operators we shall give here only 
the expressions for these quantities: 



^-k + M (t) ^i (t) = y^ 11 <w (t) v (t) - f k 2) * ( p'* )a W (t) a -i 

(2.11) 

+ (il)B a (t) - g^ 2) *(£)B_ £ (t) + 

and similar expressions for the other combinations. 
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Such an expansion (2.10) of a Heisenberg field in terms of free fields 

has been denoted * 13 ^ by the name of dynamical map. The dots 

in (2.11) stand for higher-order normal products. Let us note 

that the terms of the expansion (2.11) are linearly independent 

of each other; therefore the problem of convergence does not 

exist. In (2.10) the time dependence of a k+ + comes only 

through the time dependence of the operators a fc+ ^ and B k : 

-iE. t 

<w (t) - \ M e (2 - 12) 

-iw,t 

B k (t) = B k e k (2.13) 

The quantities 8 k , f (1) , f (2) , g U) , g (2) , E fc and u> k are to be 
determined by requiring that ijj satisfies the canonical equations 
of motion which correspond to the Hamiltonian (2.1). 
Introducing 

H g = A | d 3 x[^ i (x)^^(x)+^(x)^*^(x)] (2.14) 
where A, the order parameter, is defined as 

A = X<0|* + (x)^ (x) |0> (2.15) 
and rewriting the Hamiltonian (2.1) as 

H = (H 0 +H s ) + (H 1+ H c -H s ) = H; + H. nt (2.16) 

the coefficients in (2.10) are fixed by requiring the 
diagonalization of H^. It is well known that 

cos2 9 k = E k /E fc sin26 k = -A/E k (2.17) 
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where 



_ = _1_ ( ^2 2 



(2.18) 



r 2 J _ A 2,l/2 
E k = Ie k +A ] 



(2.19) 



Introducing the notation 



(2.20) 



we shall define the following Bethe-Salpeter amplitudes 



F(x,y;p,q) = <0|T[$(x)$ t (y) ] |a p ^a qf > (2.21) 



G(x,y;?,) = <0|T[$(x)$ + (y) ] |B Jl > 



(2.22) 



where | aa> denotes a state of two quasi-electrons : I ct q+ > ~ 
oip + ci^ | 0> and | B> a state of single boson | B^> = B^|0>. By 
requiring that the coefficients of the dynamical map have the 
form of retarded functions, the quantities f^, f^, g^ 1 * 
and g^ can be expressed as ^ : 



f k 1} (p,*> 



_ E _ ie [cosO k cos9 k _ Jl F (1) (p,il)+sin6 k sin9 k-Jl F v "' (p,£) 

(2.23) 



(2) 



E k +E k-£ 



- |(sin6 k cos0 k _ Jl +cosO k sin9 k _ jl ) f (p,£) 1 



f k 2) (p,Ji) 



-A 



[sinO k sin9 k _ £ F (1) (p , £) +cos g^osg^^F v " (p,£) 

(2.24) 



(2) 



E k +E k-£ +E+ 



+ i (sin9 k cos9 k _ Jl +cos9 k sin6 k _ jl ) f (p , fc) ] 
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(2.25) 

- 2 (sin8 k cos9 k _£+cos8 k sin9 k _ £ ) g U) ] 



^k 2) <*> = E V+ E V [3in \ SinG k-£ G(1) (^)+cose k cos9 k ^G (2) U) 

k k_£ £ (2.26) 

+ i(sin9 k cose k _ Jl +cos0 k sine k _ Jl )g(il) ] 

In (2.23) and (2.24) E = E +E . The functions F (1) , F (2) , f, 

P P *- 

G (l) ^ G (2) and g are related to the B _s reduced amplitudes 

F(x,x;p,£-p) = e i( **^~ Et) F(p,JU (2.27) 
i (x- t-<o ff t) 

C(x,x;£) = e % G(£) (2.28) 

by 

F (1) = F 12 F (2) = -F 21 f = F 22 - Fll (2.29) 

G (1) = G 12 G (2) = -G 21 g = G 22 - Gll (2.30) 

The expressions (2 . 23) - (2 . 26) are valid in the pair approximation 
in which we take into account only the Feynraan diagrams which 
conserve the number of pairs. This approximation will be used 
consistently in all of the remaining part of the present work. 
It should be noted that in this approximation the terms denoted 
by dots in (2.11) do not contribute. 

In the absence of a Coulomb potential the Bethe-Salpeter 
equations for the amplitudes (2.21) and (2.22) have been derived 
in Ref. 2. When the Coulomb interaction is taken 
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into account these equations take the form: 

D(x,y) = D Q (x,y) - i j d 4 zS (x-z) F (z) S (z-y) 

(2.31) 

where 

D(x,y) = <0|T[$(x)$ f (y) ]| l> (2.32) 

| l> being any physical state different from the vacuum, i.e. 
<0|£> = 0. D q is the B-S amplitude in the absence of inter- 
actions and S (x) satisfies the equation 

[e(V 2 )T 3 + - A Tl JS(x) = 16 (x) (2.33) 

where t ^ and x ^ are the Pauli matrices. The function F(x) in 
(2.31) is defined as 

F ll (x) = _ ^t D 22 (x,x) + \ c l(^) [D 2 2 (x,x} " D ll (x ' x) ^ 
F 12 (x) = XD 12 (x,x) F 21 (x) = AD 2i (x ' x) (2.34) 

F 22 (x) = -A{D i;L (x,x) + \ q(t) [D i;L (x,x)-D 22 (x,x) ]} 
q(l) is given by 

qi t) «^sL[ d 3 x-i-e iJ * x - e l x l = |3lsL ( 2.35) 

A J \x\ \U 2 +t 2 ) 

where the infinitesimal quantity c has been introduced in order 
to define integrations which contain Coulomb potential. From 
(2.31) it is straightforward to derive the equations for the 
B-S reduced amplitudes (2.27) and (2.28): 



F (1) (p,A) = cosB cos9 0 +Q 1:L (t,E)F (1) (p,ii)+Q 33 (t,E)F (2) (p,J! 

p p X. 



- G; <t,E) Il+q(t)Jf (p,H) 

F (2) (p,&) = sine sin9 r ^_ 5 +Q 33 (l / E)F (1) (p,£)+Q 22 (I,E)F (2) (p,i 

p p x. 

- Q 23 (t,E) [l+q(t)]f (p,l) 

(2.36) 

f (p,A) - cose sinG „ + sin9 cos9^ -2Q 13 (1 ,E) F ( 1} (p,£) 

r p p-J!, p p-x, 

- 2Q 23 (t,E)F (2) ( P ,£) + IQ 12 (1,E)+Q 33 (1,E) J [l+q($)]f (p, 
for the quasi-f ermion B-S amplitude, and 
( G (1) (1) = Q 11 (l, t o jl )G (1) (JD+Q 33 (1,oj^)G (2) (A) 

- Q 13 (l,oa a ) [l+q(l) JgU) 
G (2) (A) = Q 33 (1,03 £ )G (1) (£)+Q 22 (l,o) (i )G (2) (A) 

- Q 23 (l,^)[l+q(l)]g(£) (2.37) 

gU) = -2Q 13 (1,03 A )G (1) (£)-2Q 23 (t,w £ )G (2) (A) 

+ IQ 12 (l,o^) +Q 33 (1,(^)1 [l+q(t)]gU) 
for the boson B-S amplitude. 

The coefficients Q 1 ^ are defined in the Appendix A. It is 
interesting to note that the effect of the Coulomb potential 
can be conglomerated only as a factor q(^) multiplying the 
amplitudes f(p,&) and gU). We shall postpone till the later 
sections the problem of solving the two set of equations (2.3< 
and (2.37). Here we shall limit ourselves to mention that th< 
equations should be solved taking into account the self-consii 
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condition 

1 + * | 7777 SSf: - 0 < 2 - 38 ' 

' (2ir) k k-£ 

This equation is the generalized gap equation proposed and 
justified in Ref . 2. By means of this equation the following 
important relation can be derived: 

1 - Q + (1,E)-Q 33 (£,E) = (E 2 -^)R(t,E) (2.39) 

where Q + is defined in the Appendix A and 

E +E 

, 3 k+±£ k-±l , 

R( t /E) = x " 4E ^ E , Ri . +E - )^] (2 - 40) 

E , +e , (e , -e , ) 

-2 A_ r d k 2 2 2_ _2__ — f9 41 x 

W * = ITfTT ,o ,3 4E ' E . [(E . +E . )2- E 2] (2 ' 41) 

R(I.E) J (2TT) k+ l £ k _l £ k+ l^ ^ 



3) Calculation of the Fermion Current 
a) Charge Density 

Let us consider the charge density 

p(x) = i^(x)Tj) f (x) + \J»|(x)* + (x) (3.1) 

and let us study the contribution due to the quasi-fermions . 
From the definitions (2.21), (2.27) and (2.29) we have 

<0|p(x)|a a £ > = f(x;p,£) (3.2) 
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so that in the pair channel: 



J (2tt) j (2tt) 

(3.3) 

Solving the system of linear equations (2.36) f (p,£) is found to 
be 



f(p,A) = k [sin(8 + 0„ .)X,(t f E)-2cos(0 +6 ff )X ? (t f E) 

W($,E) P P * 3 p p x, ^ 

- 2Ecos(0 -0 a )Z 0 (t,E)J (3.4) 
where E = E +E_ „ and 

W(t,E) = {R(l,E)Y(t f E)-Z (I ,E)+q(£) [Y * (£,E)-Z' ($ ,E) ] } [E 2 -W 2 K (£ ,E) ] 

(3.5) 

The quantities that appear in (3.4) and (3.5) are defined in 
Appendix B. From the expressions (B.12), (B.13) and (B.15) 
W(1,E) takes the form 

W(£,E) = R 2 d,E) (E 2 -w 2 -4A 2 ) [E 2 - W p(l) ] (3.6) 
where to 2 (£) is defined by 

P 

co 2 (l) = a) 2 [l+4A 2 R+4A 2 Rq(t) ] (3.7) 



Substituting (3.6) in (3.4) and using the expressions (B.10) , 
(B.ll) and (B.14), we obtain the following expression for f(p,£): 
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sin(8 + 9 J (E 2 -ok)+2AEcos(0 -9 -) 
f(p,A) = P P ~* ^ j-* 



E -a> (1) 
P 



Let us now denote the quasi-f erraion field by 

* = 



V 
♦I 



which satisfies the equation of motion 



[e{V)T 3 + - A Tl ]$(x) = 0 



Then let us define the following operators 



3(0) = i[*^* f -(^)* t +*^**" ( M } *+ J 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



P (0) = *IV*I*+ 



(3.12) 



The expressions (3.11) and (3.12) correspond to the current 
and charge density in the usual Hartree-Fock approximation. 
In the self-consistent method they correspond to the current 
and charge density calculated to the lowest order approximation 
that is to the first order in the dynamical map. It is well 
known that at such order the current is not conserved, but 

a^ 0 ' = M (0) ♦ -ft p (0) = 2i4<* t *+-*X> (3 - 13) 



By diagonalizing the equation (3.11) we find the express 
of <{>^ and <J>^ in terms of free fields: 



ion 



f (x) = j^—j e lx *^[cos8 k a k1 . (t)-sin9 k at ki (t) ] (3.14) 



<j, + (x) = |^JL- T e lx *^[cose k a k4 (t)+sine k a! kf (t) ] (3.15) 



Substituting the expression (3.8) in (3.3) and introducing 
the Fourier transforms of (3.11) and (3.13) we can easily see 
that the Fourier component of the fermion contribution to the 
charge density has the form: 



> F (jfc,t) 



+a)2(t) 
P 



r Z 2 -2 



P (0) (l,t) 



By noticing that uL vanishes at 1 = 0 and that 



(3.16) 



9t p " Vv 



= 0 



£=0 



(3.17) 



we obtain from (3.16) the important result 

|d 3 x p F (x) = 0 



(3.18) 



b) Current Density 



The electron current density is defined by 



(3.19) 



Introducing the dynamical map (2.10) and taking only the 
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terms that contribute in the pair approximation we obtain 



■* t * lfd^kd^fl, ix«jt* . , rt A . - , , , - 

D(x) = - - 7- e ksm(9 -. -9 , ) a , (t)a , 

m J (27T) b k+±& k-±Jl -k+±JU k+±JW 



(t) +h.c. 

_j j. 

2 2 2 2 

(3-20) 



-t 

Next introducing the explicit expressions of a and a in terms 
of a, a" 1 " and considering the quasi-f ermion contribution in the 
pair channel we have 

f = 3<0) + Jd) (3 . 21) 
where ] (0) is defined by (3.11) and is given by 

3 <l 'w = l777Tfc elS ' %(5 ' 1)o ^p* lt) V (t)+h - c - 

J (2tt) j (2tt) 



(3.22) 

$(p.t) - if^S« 8in(B -I -6 , )If U) 1 (p,*)-f (2, 1 (p,JL)] 
m J (2tt) k+ijl k-jfc k+±& k+±£, 

(3.23) 

By means of (2.23) and (2.24) we can express 5(p,t) in terms 
of the B-S amplitudes: 

t r^3, t t Asin(9 -9 ) m 

J(P,1) = J tP 3J S ± r-^T {{E + +E - )cos{e +- 8 - ) IF <P'*> 

£ 2 J(2tt) J m <E + +E_) 2 -E^ + (324) 

+ F (2) (p,£) ]+Ecos(6 + +9_) [F (1) (p,Jt)-F (2) (p , it) ] -Esin ( 9 + + 9_) f (p , %) } 

where by obvious notation E = E , ,9 + = 9 , and E = 

k±±l ~ k±|l 

E +E „. Calculating J at first order in I and by means of 
P P-A 

(B.16) and (3.8) we have 
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J(p,t) = 2AR($,E) [F (1) (p,£)+F (2) (p,H)] ~-l (3.25) 



(2) 



-2 



where the definition (2.41) has been used. From the B-S 
equations (2.36) we have 

(1-Q + -Q 33 ) (F (1) +F (2) )-EQ-(F (1) -F (2) )+2EL"[l+qd)]f = cos(8 p -8 ) 

(3.26) 

neglecting the term proportional to f^-F^ 2 * since it is of 
higher order in 1 than the present order of calculation and 
using the explicit expression (3.8) for f(p,t) we find 



F (1 > + F (2) = 



(t) 2 



= - 0 - — £ r—cosO n -8 n 0 ) 

R[E 2 -c d) 2 \ ~ 2 P P ~ l 



UK 



P ^ 

+ 2AER[l+q(t) ]sin(9 + 9 . ) 

p P — *> 



(3.27) 



where Eq. (2.39) and definition (3.7) have been used. From 
(3.25) and (3.27) we obtain for the Fourier transform of ]^ : 



2 



. c. 



+ w (I) 2 
9t 2 P 



*4-|r P (0) (l,t)-^ V a) 2 t^ (0) (t,t> 

^2 3t P 



(3.28) 



where 3 (0) and p (0) are defined by (3.11) and (3.12), 



respectively. Substituting (3.2 8) in (3.21) we get the 
following expression for the quasi-f ermion contribution to the 
current density: 
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3 F (t,t) = 



P 



at 



3t 2 1 



-2 



(3.29) 



By noting that at the present order of computation 



(15) 



(3.30) 



where v = v //T, and that 

O F 



4A^RqU) u % 



a 2 

4Tre n „ 
2 s 2 

= w 

m p 



(3.31) 



where u p is the plasma frequency, we can rewrite the charge 
and current density in the coordinate space as 



j F (x) 



Q6 -3 3 

I 1 1 IM II 



u\> u v .(0) 
2/ 2 



(3.32) 



where the symbol Q stands for the "D* Alambertain operator" 



□ 



2 . . 2 

V dt 
O 



- V 2 



(3.33) 



with a velocity of propagation equal to v . The conservation 
of the fermion current 

3 j F (x) = 0 (3.34) 



is immediately verified from the expression (3.32). 



4) Calculation of the Boson Current 

In this section we shall compute the contribution of the 
boson field to the charge and current density. In order to 
obtain a solution for the dynamical map (2 . 10) - (2 . 11) we must 
solve the B-S equations (2.37). Let us start by noticing that 
this system of linear equations will admit a non-trivial 
solution only when 



where W is given by (3-6) and w £ is the energy of the boson 
field B. The equation (4.1) will fix the energy spectrum of 
B; from (3.6) we find 



It should be noted that these expressions are valid 

neglecting terms of the order of 0(£ 4 ). Therefore we find that 

the dynamical map admits a solution corresponding to the 

presence of two collective modes among the free fields. The 

presence of a collective mode whose energy spectrum is given 

by (4.2) has not been reported in preceding studies (Cfr. 

Refs. (1) and (3)). We have not examined in detail this solution 

and we defer to a future article the problem of studying if 

this collective mode will contribute to dynamical effects. Only 

we would like to comment that this field does not contribute 

to the generator of the gauge transformation, so that its 

presence does not affect the boson method as described in 



WU,<d 0 ) = 0 



(4.1) 




(4.2) 



(4.3) 
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Ref . (4) . In the present article we shall limit ourselves to 
consider only the boson field with an energy spectrum given by 
(4.3). According to (3.7) and (2.35) and by noticing that 
4A 2 R(0,0) = AN(0), where N(0) is the density of states at the 
Fermi level, oi^ can be expressed as 

2 2 

9 9 9 7 e % 
oj* = v £ [1+AN(0)] + or (4.4) 

* ° JT+e 



At first let us note that in the limit e ■* 0 (i.e. no 
Coulomb potential is present) we find a collective mode with a 
sound velocity 

v = v [1+AN(0)J 1/2 (4.5) 
s o 

already found by several authors (1 ' 2 ' 14 * . It is worth noticing 
that this expression agrees with the result found in Ref . (14) , 
but does not agree with the result of Anderson ^ which had a 
factor of 4 multiplying the second term in the bracket. 

In the charged case special attention must be devoted to 
the quantity e which, roughly speaking, corresponds to the 
inverse of the linear dimension of the system, e - 1/L. Keepina 
e as an extremely small but finite quantity corresponds to 
the case when the system is not really infinite. From (4.4) 
we obtain: .' 

for K<e (4.7) 



-I ■ < \ 



w 2 +v 2 & 2 for l»e (4.6) 



This form for the energy spectrum of the boson field, already 



reported in Ref. 3, is in contrast to previous s tudies ^ * 
which found that in the presence of the Coulomb potential the 
collective mode acquires a finite mass which coincides with 
the plasma value. From (4.4) we find that practically the 
energy of the collective mode starts at the plasma frequency, 
but for the presence of (4.7) it is possible to have condensation 
of bosons in the ground state. This result furnishes a proof 
of the fact that the Goldstone theorem may be true even in the 
presence of long range interactions. In section 5 we shall 
discuss the Goldstone theorem and the rearrangement of the 
symmetry in the light of the results derived. 

Let us now introduce the explicit expression for the boson 
field in terms of annihilation and creation operators: 

-it. 



B(x) = Z 



A- 



il'X-ioipt . 



♦x+iui^t 



(4.8) 



This expression reduces to the usual expression 
free field when e -* 0. The field B will satisfy the equation 
of motion 



j 2 2 

3^ 2*2 2 G % 

— r - v V +o) t~Tt- 

2 s p 2 -^2 



3t 



e -V 



B(x,t) = 0 



(4.9) 



which can be rewritten in the form: 



2*2 



3t 



5 - v^v 
s 



B(x,t) 



4tt 



d 3 y 



- £ !*-^ 2 B(y,t) 



|x-y| 



= 0 
(4.10) 



The last term in (4.10) represents the contribution due to the 
Coulomb effect. We shall also introduce the canonical conjugate 
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of B, ir: 
tt(x) = -iZ 



v I 
s 



/2 



0), 



iX-x-ioj t , -it' 
B A e *■ - Bj e 

(3) 



x+io)^t 



It should be noted' that tt is related to B by: 



(4.11) 



tt = B in the neutral case 



(4.12) 



1T(x,t) = 



1 f 3 e " y l x " y l ->2- -> 

Ttv J d y T^-f | V B ^ Y,t ^ in the char ?f ed case 



(4.13) 



where u = u /v . 

p s 

In order to solve the B-S equations (2.32) we shall fix 
the normalization of the amplitude (2.28) by requiring ^ that 
the field B satisfies canonical commutation relations. We thus 
find the normalization condition: 



(2tt} J k 



2 - |g< 2) U) | 2 } = i 



(4.14) 



where the functions g (1) and g (2) are defined by (2.25) and 
(2.26), respectively. 



a) Calculation of the Charge Density. 

From the definitions (2.22), (2.28) and (2.30) we have 

<0|p(x) |B ? > = g(x;£) (4.15) 

so that the boson contribution to the charge density is given by 

fd 3 £ ix-t-iw t 
= J e * gU)B 0 +h.c. (4.16) 



p (x) 



(2tt) 



Solving the equations (2.37) with the condition (4.14) we 
find 



g(i) = -12Av sJ y|^- (4.17) 

Substituting this result in (4.16) and taking into account the 
definition (4.11) we obtain 

p B (5,t) = -n(3)ir(x,t) (4.18) 
where the coefficient n (3) is defined by 

n(3) = -2A^R(9) (4.19) 

b) Calculation of the Current Density. 

Substituting the dynamical map (2 . 10) - (2 . 11) in (3.19) 
and considering the boson contribution in the pair approximation: 



3 Cx,t) = - 



fd 3 k (d\ is^-v, 

7773 77T3 e ksin(9 , .i„~ 9 , i. 

(2tt) j 



(2tt) J k+±S, k-±& 

+ [g (1 ! U)-g (2 ! (£)]B +h.c. (4.20) 
k+^S, k+±£ 

From (4.17) and from the B-S equations we can evaluate the 
functions g^ and g^. The final result is 

3 B (x f t) = v 2 n(3) VB(5,t) (4.21) 

where use has been made of definition (4.8). By noticing that 
from (4.9) the boson field satisfies the equation 



3 2 2±2 2 

- V V +W 

_9t 2 3 P 

we obtain from this relation and from (4.13) the important 
result: 

*(x,t) - v 2 V 2 B(x,t) (4.23) 

Owing to this relation the conservation of the current is 
immediately verified: 

Vy = ^ B + "It pB = n(3) I v s^ 2b ^ 3 = 0 (4 ' 24) 

As we have already mentioned Eq. (4.22) is very 
different from the equivalent equations obtained in other 
approaches. The important point in (4.22) is that this equation 
admits the solution 

$ 2 B = 0 (4.25) 

The existence of this solution will permit boson condensation 
in the ground state which recovers the gauge invariance of the 
theory. 

5) Long Range Correlations, Goldstone theorem and Rearrangement 
of Symmetries in Superconductivity . 

The conditions of applicability of the Goldstone theorem 
can be schematized in the following three points: 

a) a Lagrangian theory which has an internal symmetry 
group ; 

b) asymmetric nature of the vacuum (i.e. conditions of 
long range order) ; 



V B(x,t) 



(4.22) 



c) the generator of the symmetry group is time 
independent^ 6 ^ „ 

The last condition has brought a series of interesting 
speculations concerning the study of the cases in which 
the Goldstone theorem appears to be invalidated. In particular, 
since the condition c) is always satisfied in the case of 
manifestly covariant theories, attention has been devoted to 
non relativistic theories and to relativistic theories which 
do not present manifest covariance. The experimental evidence 
that long range excitations are present in superf luid liquids , 
but absent in superconducting systems, and the theoretical work 
of Anderson brought people to speculate about the possibility 
that the Goldstone theorem breaks down when long range fields 
are present in the theory. This circumstance has special 
interest in the physics of elementary particles if one wants to 
explain the presence of massive observed particles without 
appealing to the introduction of breaking Lerms in the 
Lagrangian. 

The common argument that relates the presence of long 
range fields to the invalidation of the Goldstone theorem is 
the failure of the condition c) . If i is the conserved current 
we have from the conservation law 



The generator N will be time independent if and only if the 
surface contribution of the current vanishes. In non relativist 




(5.1) 



theories the presence of long range forces may imply that 

there exists a flow of current passing through the surface which 

contains the system; as a consequence in (5.1) we shall have a 

contribution from the surface integral and the generator will 

be time dependent. In relativistic theories which are not 

manifestly covariant the requirement of causality does not 

necessarily apply, so that one might expect a non vanishing 

(19) 

contribution from the spatial integral in (5.1) 

We would like to observe that this argument does not bring 
much information since no general statement can be reached 
about the effects of the original invariance: the problem must 
be solved case by case according to the dynamics of the particular 
system under study. Also this scheme of thinking may bring 
serious difficulties in quantum field theory since the existence 
of incoming fields is not necessarily justified in the case in 
which effects from infinity occur. As we have seen our results 
contain surface effects only through e which fixes the shape of 
the energy spectrum around I = 0; we will show that the generator 
is time independent, and does not depend on these effects. This 
circumstance preserves the original invariance of the theory. 

In relativistic theories the effects of long range fields 
on the Goldstone theorem have been studied in different models. 
The Goldstone model * 20 * interacting with a gauge field has been 
considered by Higgs in a semi-classical treatment, and by 
Englert and Brout (11) in the lowest order perturbation theory. 
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(21) 

The Nambu and Jona-Lasinio model with an axial vector gauge 

field coupled to the helicity current has been studied by 

f 12) 

Freundlich and Lurie . In all the models considered the 
results turn out to coincide with those reached by Anderson 
in superconductivity. The introduction of the gauge field 
eliminates completely the Goldstone particle (present in the 
original Golstone and Nambu and Jona-Lasinio models) while 
at the same time the gauge field acquires a finite mass. 

All these studies in relativistic and non relativistic 
theories seem to show that the Goldstone theorem is not valid 
anymore when long range interactions among the particles are 
present at the level of the Heisenberg fields. However, a more 
careful analysis shows that all these considerations are not 
completely satisfactory and many problems remain unanswered. 
We have already mentioned that the original invariance of the 
theory must be conserved even at the stage of the physical 
particles; spontaneous breakdown of symmetry implies in reality 
that the form of the invariance transformation takes a different 
shape when written in terms of the physical fields. If we examine 
all the works previously mentioned it is easy to see that the 
original invariance transformations vanish entirely from the 
theory when expressed in terms of physical fields; in other 
words the physical fields do not contribute to the generator 
and remain completely unchanged under the transformation. This 
conclusion is rather surprising since we should expect that at 
the level of the physical fields some entity must carry the trans- 
formation. This aspect is closely related to the fact that apparently 
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the Goldstone particle has completely disappeared from the theory 

without leaving any trace . When the rearrangement of the 

symmetry has been studied in different models using the self- 

(13) 

consistent method/ it has always been found that at the level 
of the physical fields the noldstone particles are the only fields 
which contribute to the generator and which transport all the 
symmetry attributes. 

Summarizing, we would like to say that the problem of 
spontaneous breakdown of symmetries when long range correlations 
are present has not been studied with sufficient care, and more 
analysis is required. One of the motivations of the present 
work has been to study this problem in non relativistic theories 
taking in consideration the BCS model in superconductivity. By 
using the self-consistent method we have solved the model in 
the pair approximation and an answer can be given to the questions 
asked before. 

Let us consider the constant gauge transformation of the 
electron field 

-iNO, iNO i9. , c 
^ y ifj' = e ij>e = e i|j (5.2) 



where 



N = Jd 3 xp(x) ; (5.3) 



0 is a constant real parameter and p is defined by (2.6). The 
transformation (5.2) leaves the Ilamiltonian (2.1) invariant. 
Let us now examine how the physical fields J> + ^ and B change 
under the transformation induced by (5.3) . Taking into account 
the results (3.18) and (4.18) the generator (5.3) can be 



expressed in terms of physical fields and takes the form 

N - -n(3) fd 3 XTT(x} (5.4) 



The first consequence of this expression is that the quasi- 
fermion fields ^ do not contribute to the generator and 
remains invariant under the gauge transformation: 

* V* (5 - 5) 

while the boson field transforms as 

B(x) — ► B'(x) = e ~ iN9 B(x)e lN9 - B(x) + n9 (5.6) 

IT *■ TT (5.7) 

where r\ = r\(Q) . The second important consequence which follows 
from (5.4) is that the generator is time independent 

N = 0 (5.8) 

This result is immediately proved by recalling the expressions 
(4.11) and (4.4). The conclusion (5.8) can be realized without 
introducing Goldstone excited levels by ueans of the presence 
of the Laplace operator in Eq. (4.22); owing to the presence 
of V 2 the system exhibits a condensation in the ground state 
and the transformation (5.6) leaves invariant the equations of 
motion (4.9) and (4.10) for the boson field. We thus find that 
the gauge transformation (5.2) of the electron field is re- 
arranged at the level of the free fields in the shape of a linear 
transformation of the boson field; the original gauge invariance 
of the theory is not lost but is recovered by the presence of 



the boson field which acts as a Goldstone particle. From (5.5)- 
(5.7) it also follows that the electron fields ^, expressed 
in terms of physical fields, take the form 

ii B 

^ + = e n F(<j>,VB,Tr} (5.9) 

where F is a certain unknown function. 

In conclusion, we would like to state that our results, 

f 1) 

in contrast with previous studies , seem to indicate that the 
Goldstone theorem is still valid in the presence of long range 
forces in the sense that a boson condensation is still permitted. 
A collective mode whose energy goes to zero in the limit of 
of long wave lengths is present in the theory? this massless 
particle is not only responsible for the asymmetry of the 
ground state, but rearranges and preserves the symmetry of the 
theory when formulated in terms of quasi-f ields . 

(22) 

Similar conclusions have been recently reached in 

(23 

relativistic theories in the study of non convariant solutions 

(11) 

of the generalized Nambu and Jona-Lasinio model 
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Appendix A 



ii 

The kernels Q J which appear in the Bethe-Salpeter 
equations are defined by: 



Q 11 = Q + + EQ~ Q 22 = Q + - EQ" 



J C2tt ) 
3 

Q"(£,E) - -|f^"^( E k E k Hl +e kHlV f <M,E) 



Q 13 = EL" + L + Q 23 = EL" - L f 



L + (1;E) M d k 3 (£ k +e k _ £ ) (E k +E R _^) f (k ,1 ,E) 
J (2tt) 

L~(£,E) = -AR(£,E) Q 33 (t,E) = -2A 2 R(1,E) 

Q 12 (1,E) = -IIt^T^E^-^^) (E k+ E k _ A ) f (S.t,E) 
J (2tt) 

where R(£,E) is defined by (2.36) and 
f (k,£,E) 



F 'k E k~£ [(E^E^) 2 ^ 2 ] 



Appendix B 



The quantities K, X 2 , X 3 , Y, Y r , Z, Z', Z 2 introduced in 
section 3 are defined in terms of the coefficients Q 1 - 1 in 
Appendix A. Their explicit expressions are 



K= [RY+q (H) Y' ] [RY-Z+q (t) (Y ' -Z 1 ) ] ~ (B.l) 

X 2 = L + (1-Q + -Q 33 )+E 2 L"Q" (B.2) 

X 3 = (1-Q + -Q 33 ) (1-Q + +Q 33 )-E 2 Q" 2 (B.3) 

Y= (1-Q + +Q 33 ) (1~Q 12 -Q 33 )-4L +2 (B.4) 

Y * = -R(1~Q + +Q 33 ) (Q 12 +Q 33 )-4RL +2 (B.5) 

1 ± = Q""(1-Q 12 -Q 33 )+4L + L~ (B.6) 

Z 2 = L~(l-Q + +Q 33 )+L + Q" (B.7) 

Z - Q - Z 1 +4L~Z 2 (B.8) 

Z* = Z-Q" 2 (B.9) 



For values of I sufficiently small with respect to k , the 
Debye momentum, the domain of integration in Eq. (2.34) can be 
considered 5,-independent ; using the relation (2.35) and the 
expressions of the coefficients Q 3 the following results can 
be obtained: 

X 0 (t,J2) * OU 4 ) (B.10) 
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X 3 (i,E) = R 2 (E 2 -cjJ) (E 2 -(dJ-4A 2 )4-0(£ 4 ) (B.ll) 
Y(1,E) = R(1+4A 2 R) (E 2 -w 2 -4A 2 )+0 (A 4 ) (B.12) 



2($,E) = Z'(~£,E) = Y'($,E) = 4A 2 R J (E 2 -wJ-4A' ! )+0(jl 4 ) (B.13) 

Z 2 (£,E) = -AR 2 (E 2 -5 2 -4A 2 )+0(£ 4 ) (B.14) 



Substituting these values in the expression of K we obtain 

K($,E) = l+4A 2 R(t,E)+4A 2 R(i / E)q(t)+0 (I 2 ) (B.15) 
The result (B.10) is derived observing that 



AQ--L+ - f- 



fd\ £-1 E k E k-& „ ft „4. 



(B.16) 
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